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Abstract 

We discuss DGLAP and BFKL evolution equations in the N = 4 supersymmetric 
gauge theory in the leading and next-to-leading approximations. Eigenvalues of the BFKL 
kernel in this model turn out to be analytic functions of the conformal spin, ft allows us 
to find the residues of the anomalous dimensions of the twist-2 operators in the points 
j = 1,0,—1,.... from the BFKL equation in an agreement with their direct calculation 
from the DGLAP equation. The holomorphic separability of the BFKL kernel and the 
integrability of the DGLAP dynamics in this model are also discussed. 

PACS: 12.38.Bx 


1 Introduction 


The Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation 10 is used now together with 
the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [[| for a theoretical 
description of structure functions of the deep-inelastic ep scattering at small values of the 
Bjorken variable x. In this kincmatical region the structure functions are measured by 
the HI and ZEUS collaborations [|]. The radiative corrections to the splitting kernels for 
the DGLAP equation are well known |J. Although the BFKL equation in the leading 
logarithmic approximation (LLA) was constructed many years ago, the calculation of 
the next-to-leading corrections to its kernel was started only in 1989 |6] and completed 
comparatively recently 0- 01- 

In supersymmetric gauge theories the structure of BFKL and DGLAP equations is sim¬ 
plified significantly. In the case of the extended iV = 4 SUSY model the next-to-leading 
order (NLO) corrections to BFKL equation were calculated in ref. 0 for arbitrary val¬ 
ues of the conformal spin n in a framework of the dimensional regularization (DREG) 
scheme. In the section 3 below the results are presented in the scheme of the dimensional 
reduction (DRED) which does not violate supersymmetry. The analyticity of the eigen¬ 
value of the BFKL kernel over the conformal spin n gives a possibility to relate DGLAP 
and BFKL equations in this model, as we show below. Further, the eigenvalue in both 
schemes (DREG and DRED) can be written as a sum of the eigenvalues of holomorphic 
and anti-holomorphic operators (see Section 3). 

Let us to introduce the unintegrated parton distributions (UnPD) (p a (x , k 2 L ) (hereafter 
a = q,g,p for the quark and gluon densities respectively) and the (integrated) parton 
distributions (PD) f a (x,Q 2 ), where 

fa(x,Q 2 ) = f dk 2 ± (p a (x,k 2 ± ). (1) 

Jk 2 ± <Q 2 

The DGLAP equation relates the (integrated) parton distributions having different 
values of Q 2 . It has the form: 

TT7vM x ’Q 2 ) = -1U/«M 2 )+ ^(x/y)My,Q*), 

d In Q z J x y ^ 

W a = Y,! 1 dxxW a ^{x\ (2) 

b Jo 

where the last term in the r.h.s. is the Mellin convolution of the inclusive transition 
probabilities W b _> a (x) and PD f b (x,Q 2 ). Usually the first and second terms in the right- 
hand side of the equation are unified by modifying W b -> a in the form of the splitting kernel 
W b ^ a . 

■ i d n2 fa(x,Q 2 ) = I —J2W b ^ a (x/y)f b (y,Q 2 ) (3) 

dmQ z Jx y 

It is known, that DGLAP equation (3|) is simplified essentially after the Mellin trans¬ 
formation to the f-channel angular momentum j representation: 

rf li^g 2 /a(j’Q 2 ) = 5^7 ab(j)fb(j,Q 2 ), (4) 
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where 


fa(j,Q 2 )=[ dx x J 1 f a (x,Q 2 ) (5) 

Jo 

are the Mellin moments of parton distributions. The Mcllin moment of the splitting kernel 

7 ab(j) = I dx X 3 ^' 1 W b ^ a (x) (6) 

Jo 

coincides with the anomalous dimension matrix for twist-2 operators (k These operators 
are constructed as bilinear combinations of derivations of the fields describing the partons 
a (see Eq. (|D below). 

The BFKL equation relates the unintegrated parton distributions having various values 
of the Bjorken variable x and has the form: 

dhitl/x) tP9 ^ X,k ^ = 2 o} ^~ k ^ (p si x ^ k ±) + j d 2 k , ± K(kx,k' ± )tp g (x,kl), (7) 

where uj(—k\) is the gluon Regge trajectory. 

Let us introduce the local twist-two operators: 


O 9 

Ml,- 


_ C F) D 

PP1 ^2 ^P-2,' 

D C 1 

" ^pj-i^ppji 

d 9 

Ml,- 


= SG Pfll D p2 D p3 . 


o q 

Mi,- 


= S^ tll D P2 ...D j 


o q 

Ml,- 


= 5^75 7 M1 D P2 ... 


0 V 

Mi,- 


= S$D P1 D P2 ... L 



where the last operator is constructed from derivations of the scalar field <f> appearing 
in supersymmetric models. The symbol S implies a symmetrization of the tensor in the 
Lorenz indices /r 1; ..., /jj and a subsequent subtraction of its traces. 

The matrix elements of Of, ,, and Of ,, are related to the moments of parton 

Plv) Pj 

distributions in a hadron h in the following way 

J q dxx 3 ~ 1 f a (x,Q 2 ) = < h\n^...h^ 0^ l itl .\h>, a=(q,g,<p), 

J o dxx 3 ~ 1 Af a (x,Q 2 ) = <h\n^...h^ d^ i}tl .\h>, a=(q,g), (9) 

where the vector n 9- is light-like: n 2 = 0. Note, that in the deep-inelastic ep scattering 
h 9 = q + xp. 

The conformal spin n and the quantity 1 + aexpressed in terms of the eigenvalue 
oj of the BFKL kernel, coincide respectively with the total numbers of transverse and 

1 As in Ref. 0J, the anomalous dimensions differ from ones used usually in DIS by a factor (—2), 

i-e. 7obO') = (-l/2)7ffc /S (j) 
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longitudinal indices of the Lorentz tensor with the rank j — 1 + u + n. Namely, we can 
introduce the following projectors of this tensor 




/O a JC Tl J a \n\ 

/il | n | (Ti ,...,a'|n| -L -L 5 


where the complex transverse vector l± is given below 

11 = + = 0. 


( 10 ) 


Anomalous dimension matrices 7 a b(j) and lab{j) for the twist-2 operators and 

do not depend on the different projections of the tensors due to the Lorentz 
invariance. 

In the gluon case the usual light-cone projections entering in the DGLAP equation 
(see Eqs.fllOl)) are 

h^...n^ <P\Ol tl ^.\P> = I dxx J ~\f g (x, Q 2 ). (11) 

J 0 

The mixed projections 

n M ...n w+ "Z^ 1+ "...Z^' < P\0^’ ’^\P > ~ £ dxx“ J d 2 k ± (J^j (p g (x,k 2 ± ), (12) 


can be expressed in terms of the solution of the BFKL equation. Generally the corre¬ 
sponding operators have higher twists. 


Thus, it looks possible to obtain some additional information about the parton x- 
distributions satisfying the DGLAP equation from the analogous fc^-distributions satis¬ 
fied the BFKL equation. Moreover, in the extended A^ = 4 SUSY the /3-function equals 
zero and therefore the 4-dimensional conformal invariance could allow to relate the Regge 
and Bjorken asymptotics of scattering amplitudes. 


Our presentation is organized as follows. I 11 Section 2 we discuss the relation between 
DGLAP and BFKL equations in the leading logarithmic approximation. In Section 3 we 
review shortly the results of Ref. [|SJ and rewrite them in the framework of the DRED 
scheme. Section 4 contains the information about the anomalous dimensions calculated 
independently with the use of the renormalization group. A summary is given in Conclu¬ 
sion. 


2 Anomalous dimensions of twist-2 operators and their 
singularities 

In the leading logarithmic approximation (LLA) for the BFKL equation the contribu¬ 
tion of fermions is not essential and therefore in this approximation the integral kernel is 
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the same for all supersymmetric gauge theories. In the impact parameter representation 
due to the conformal invariance the solution of the homogeneous BFKL equation has the 
form (see fldfl ) 


E v ,n{]hOi P20) =< 0 (dl) Om,m(po) 0(^2) > 


P 12 P 12 

P 10 P 20 / \p* 10 p *20 


m 




(13) 


where 


m = 1/2 + iv + n/2, m = 1/2 + iv — n/2. 


are the conformal weights. Here we introduced the complex variables p k = x k + iy k and 
denoted p k i = p k ~ Pi 

For the principal series of the unitary representations the quantities v and n are 
correspondingly real and integer numbers. The projection n of the conformal spin \n\ can 
be positive or negative, but the eigenvalue of the BFKL equation in LLA 



depends only on |n[. We shall imply below that n is positive, i.e. 


n = \n 


(14) 


(15) 


Note, that eq. m has the property of the holomorphic separability and corresponds to 
the pair hamiltonian of the integrable Heisenberg spin model (see [11|-[16|). 


The solution of the inhomogeneous BFKL equation in the LLA approximation can be 
written as the four-point function of a two-dimensional held theory 


/ oo r 

duC(u, n) / d 2 po 
-OO J 


E v , n {A P2o)E/, t7l (pro : pvo) 
to — u>°(n, u) 


,( 16 ) 


where C(is,n) is expressed through the inhomogeneous term of the equation with the use 
of the completeness condition for E VyU (see [|Rj]). 

For pu — > p 2 > we obtain from the integration region Po ~^ Pi' : 

< 0(Pi)0(d2)0(Pu)0(p?) > ~ 


where is a solution of the equation 

u = a;°(n, v) (18) 


E 


dvCW, n) 


U — cu°(?z, v) 


\ Ev u ,n(PlV 1 P21') I 
2 ^C (u w ,n) -- 7 -Pi' 2 ' 


1+2 iv u t Pl'2' 

\Pl'2'. 


n/2 


(17) 
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with Imvu < 0 . 

The above asymptotics has a simple interpretation in terms of the Wilson operator- 
product expansion 


lim 0 (pn) 0 (pir) 

Pi' > P 2 ' 


E 


Cjy^^n) 
a y(n, v w ) 



P 1 ' 2 ' 
PV 2 ' 



(19) 


where 

r w = - + iVui (20) 

is the transverse dimension of the operator 0 Vui , n (pii), calculated in units of squared mass. 
This operator is the following projection 


Ov u ,,n(Pl 1 ') 


n w ...n M+ " I?-1? O^, 




( 21 ) 


of the gauge-invariant tensor with 1 + u + n indices. 

The anomalous dimension 7 ( 0 ;) obtained from the BFKL equation in LLA has the 
poles 

Ini . , . .. q 2 N r 

r„ = l + y-i7(A 7(j)Uo = !^. ( 22 ) 

The canonical contribution 1 + |n|/2 to the transverse dimension T w corresponds to 
the local operator 

G PI j. 1 j ] 2 • •' Dl L D n ■ ■ - D a n G'pmi+u, > ( 23 ) 

because in the light-cone gauge = 0 the tensor 

GpnGpn+^n^n^ 2 ~ Ajdp 2 A^ h Ml h M2 (24) 

has the transverse dimension equal to 1 . 

The above local operator 0 Vw , n for |n| > 0 should have the twist higher than 2 because 
its anomalous dimension is singular at uj — > 0. Indeed, such singularities are impossible 
for the twist -2 operator, because for n > 0 and small ui the total number of its indices 
is integer and physical. The only way to obtain some information about the twist-2 
operators is to continue analytically the anomalous dimension 7 ( 0 ;) to negative integer 
points 


n 


— r — 1 , 


(25) 


where r is a positive integer. Then in the limit cu —» 0 the total number of indices 
j — 1 + ou + \n\ tends to — r and therefore the pole of 7 ( 0 ;) in c 0 can be interpreted as the 
pole 1 /(j + r) for a non-physical value of j for the twist -2 operator. 

In LLA one can obtain after the analytic continuation of cu°(|n|, v) 


7 (j)E-r 


g 2 N c 1 

47 t 2 j + r 


(26) 
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We remind, that for the BFKL equation in LLA the fermions are not important, but 
generally they give non-vanishing contributions to the residues of the poles of 7 (j) in 
the DGLAP equation even in LLA. Therefore the above result for 7 (j) can be valid only 
for a definite generalization of QCD. It is known [|9[, that only for the extended N=4 
supersymmetric Yang-Mills theory the anomalous dimension, calculated in the next-to- 
leading approximation for the BFKL equation, can be analytically continued to \n\ = 
—r — 1. Therefore it is natural to expect, that the above result for 7 ( 7 ) in LLA is valid 
for the N=4 case. 

Indeed, using the conservation of the stress tensor in this theory to fix the substruction 
constant in the expansion over the poles at j = —r, we obtain : 


i(i) lLA = ^r(*(i) - - 1 )) 


(27) 


in an agreement with the direct calculation of 7 (j) in this theory (see 

4 ). 


1|, |19| and Section 


In the next-to-leading order (NLO) approximation it is needed to modify the above 
procedure of the derivation of 7 (j) from the BFKL equation, taking into account the 
possibility of the appearance of double-logarithmic terms leading to triple poles at j = —k. 
We shall return to this problem in our future publication. 


3 NLO corrections to BFKL kernel in N = 4 SUSY 

To begin with, we review shortly the results of Ref. [9[], where the NLO corrections to 
the BFKL integral kernel at t — 0 were calculated in the case of QCD and supersymmetric 
gauge theories. We discuss only the case of the N — 4 supersymmetric gauge theory and 
write the formulae important for our analysis. 


3.1 The set of eigenvalues 

The set of eigenvalues for eigenfunctions of the homogeneous BFKL equation for the 
N = 4 supersymmetric theory 


uj = 4 a 


X(n, 7) + ( gX(ra,7) + S {n, 7 ) | a 


has been found in || in the following form 


xUt) 
8 ( n , 7) 


2 *( l )-«( 7 + ^)-«( l -7 + |) 

- 2 4 >(n, 7) + 2 $(n, 1 - 7) + 2((2 )x(n, 7) 
6 C ( 3 ) - *"( 7 + '-) - f 1 -7 + 1 ) . 


(28) 

(29) 

(30) 
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where T( 7 ), ^>'{z) and ^"{z) are respectively the Euler T -function and its derivatives; 
a = g 2 N c / ( 167 T 2 ) is coupling constant in DREG scheme. The function <h(n, 7 ) is given 
belowQ 

rl r l-l+n/2 I" 1 / „ I 1 \ 

$(n, 7 ) = ~ J Q dx 1 + x - ( v h / ( - ) — C( 2 )J + Li 2 (—x) + Li 2 (x) 


00 r-x) k 


+ ln(x) (( 77 , + 1) — 4/(1) + ln(l + x) + 

k =1 k + n 


00 x k 


E 


k =i ( k + n ) 2 

(-i) fe+i 


(1 - (-i) fc ) 


S + 7 + ra /2 


k + 7 + n /2 


'h'(A) + 1) - T'(fc + n + 1) + (-1 ) k+1 (f3'(k + n + 1) + f3'(k + 1))) 
ty(k + n + 1 ) - T(fc + 1 )^ 


(31) 


and 


= j 


r(^) - *'(§) 

Note, that the term 


, P H (z) = - 




p\C«,7) 


(32) 


appears as a result of the use of the DREG scheme (see (29) in ||). It is well known, that 
DREG violates SUSY. The proper procedure, which is very close to DREG and satisfies 
the supersymmetry requirements, is the DRED scheme. The results in the framework of 
the DRED scheme [] can be found from (^) by the redefinition of the coupling constant 


a —> a — a + —a j 


(33) 


which eliminates the above term (^) in (^). For the new coupling constant a the above 
expression for 00 can be written in the following form 


u = 4 a 


x( n , 7 ) + 5(n, 7 ) a 


(34) 


3.2 Generalized holomorphical separability 

Following the method of Refs. [R| [L9| and using the results of the previous section 
we can present the above NLO corrections to the BFKL equation in the form having a 
(generalized) holomorphical separability. We can split the function <f>(n, 7 ) in two parts 

2 Notice that the representation of $( 71 , 7 ) contains a misprint in ||: the factor (—l ) fe+1 in the last 
sum of ( pT| ) was substituted by (—l) fc 

3 In the calculations given in |J the result obtained in the DRED scheme with 6 scalars and pseu¬ 
doscalars is equal to that for DREG with 6 + 2e scalars and pseudoscalars, where e = (4 — D)/2 and D 
is the space-time dimension. 
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$(n, 7 ) = $i(n, 7 ) + <f> 2 (n, 7 ), 


where 


and 


$ i(«,7)= E 


(j3'(k + n + 1 ) — (—1 ) k ^'{k + n + 1)) 


fc =o k + 7 + n/2 

(—l) fc (T(n + fc + 1 ) — 'h(l)) 


E 

k =0 


(k + 7 + n/ 2) 2 


$2 (n, 7 ) = ^ 
k =0 


(/3'(* + 1 ) + (-i)^'(^ + 1 )) 

k + 7 + n /2 


” (-l) fc (T(fc + l)-T(l)) 

is (^ + 7 + E 2 ) 2 


$ 2(7 + n/2) . 


Here <h 2 (n, 7 ) depends only on m = r y + n/2 and therefore the corresponding contributions 
to u have the property of the holomorphic separability. Further, for <&i(n, 7 ) we obtain 
the simple relation 


<f>i(n, 7 ) + <l>i(n, 1 — 7 ) = / 3'(7 + n/ 2 ) [T(l) - T (1 - 7 + n/ 2 )] 


+/3 / (l — 7 + n/2) [tt(l)-tt( 7 + n/2)] , 

which can be verified by calculating the residues at 7 = —r — n/2 and 7 = l+r + n/2 
and the asymptotic behavior at 7 —> 00. Therefore one has 

( 72 77 / \ 77 - 

+ g) + “ 7 + 2>J + ^ 2(7 + 2^ _ 


/% + f) 


T(l) 


T/ n \ 

^( 7 + _) 


+ $ 2 (l-7 + f)-W-7 + f) 


T(l) - T (1 


n, 

7 + ^ 


where x(n, 7 ) is given by Eq.([29|). 

Thus, we can rewrite the NLO corrections S(n, 7 ) in the generalized holomorphic 
separable form (providing that u>q is substituted by u, which is valid with our accuracy): 


<50,7) = <t>(l+ |) +')>( 1 -7 + |) - ^(p(t + |) +p(i-7 + |)) 
^0 = 4 o( 2 *(l) - *(7 + ))) - l(l “7 + ^)) 

->(7) = /? (,) (7) + ic(2) 

0(T) = 3C(3) + - 24» 2 ( t ) + 2/3< 1 >( t )(21'(1) - *(t)) • 
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(35) 

(36) 

(37) 

(38) 










oc 


3.3 The asymptotics of “cross-sections” at s 

As an example we consider the cross-sections for the inclusive production of two pairs 
of particles with mass m in the polarised 77 collision (see 0 . 0 ): 


a(s) 




1 


where a em is the electromagnetic coupling constant, the coefficient <To(s) is proportional 
to the cross-section for the scattering of unpolarized photons and cr 2 (s) describes the spin 
correlation depending on the relative azimutal angle d between the polarization vectors 
of colliding photons. 

The asymptotic behavior of the cross-sections 07 (s) (k = 0, 2) at s —> 00 leads in the 
t-channel to unmoving singularities f u (t) ~ (u — U 4 .) 1 / 2 , where 




O’o(s) 


a 2 (s) 


4a 


X(k,^)+a5(k,^) 


= 4 ax(k, 


1 — ac(k, -) 


9vr 5 / 2 


n^O 


T/2 ' (l + C(a)), 


32y 7C(3) (l„( s /»„))- 

7 r 5 / 2 / \ 

- 7 -t /2 ■ (l + 0 (a)). 

9 • 32^/77(3) - 8 (ln(s/s 0 )) 7 V ' 



Hn, t ) \ 

X(n, i) J ’ 


(39) 

(40) 


Here the symbol 0(a) denotes unknown next-to-leading corrections to the impact factors. 

Using our results (26), (27), we obtain in the N = 4 case the following values for 
X(k, |) and c(k, ±) (k = 0 , 2 ): 


x(0,1) 

= 4 In 2, 

70.1) 

= 2C(2) + 

72 . i) 

= 2C(2) + 


= 6.0348, 


y(2,-) = 4(ln 2 — 1), 


2 In 2 


11C(3)-32Ls 3 (|)-^C(2) 


16 


= 9.5812, 


2 (In 2 - 1) 


11C(3) + 32Ls 3 (^) +147rC(2)-321n2 


(41) 

(42) 

(43) 


where (see [|(| |2T[ ) 


Lsq(x) = - 


In 2 


2 -(f) 


dy. 


Note, that the function Lss(a:) appears also in calculations of some massive diagrams (see, 
for example, the recent papers [^] and references therein). 

The LO results (£11]) coincide with ones obtained in ref. [Q . As it was shown in 01 - in 
the framework of QCD the NLO correction c® CD ( 0,1/2) is large and leads to a quite strong 
reduction of the value of the Pomeron intercept (see recent analyses |23|-[25| of various 


effective resummations of the large NLO terms). Contrary to c® CD ( 0,1/2), the correction 
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c(0,1/2) is not large (c < 3 C 7 jC> (0, l/2)/c(0,1/2) ~ 2.7), which seems to support the results 
of Ref |23], where a quite large value of the non-conformal contribution to c® CD ( 0,1/2) 
was found in the physical renormalization schemes. The values of c® CD (2, 1/2) (see 0, 0) 
and c(2,1/2) are small and do not change significantly the small LO value [|2| of the 
angle-dependent contribution. 


3.4 Non-symmetric choice of the energy normalization 

Analogously to refs. 01 one can calculate the eigenvalues of the kernel in the case of a 
non-symmetric choice of the energy normalization parameter So in eq.(15). For the scale 
So = q 2 , which is natural for the deep-inelastic scattering process, we obtain 
in DREG-scheme 


u = 4 a 
and in DRED-scheme 


X(n, 7) + (jjXK 7) + S{n, 7) - 2y(n, 7 )y'(n, 7 )) a 


oj = 4a 


X(n, 7 ) + ( 5(n, 7 ) - 2x(n, 7)x'(n» 7) ) a 


where 


X '( n , j ) = -^ X ( n , 7) = —^'(7 + + T'(l-7 + ^ 


(44) 

(45) 


3.5 The limit 7 —> 0 for n = 0 


By considering the limit 7 

i) 


0 in the Eq.(l45|) we have for n 


0 (see also the analysis in 


where 


x(0, 7 ) = - + 0(7 2 ), 

7 

1 jdDREG 

nX(«,7) + £(0,7) - 2 x(0,7)x'(0,7) = - + C + 0('y 2 ), 

3 7 

jdDRED 

^(0,7) - 2y(0,7)x / (0,7) = -+ C + 0(7 2 ), 

7 


(46) 


b dreg = b dred = 0> and c = 2((3). 
3 


(47) 


Analogous to ref. 00 with the use of eqs.([4B|) and (El one can obtain the expression 
for anomalous dimensions of twist -2 operators 7 at oj —» 0 (i.e. near j = 1) 
in DREG-scheme 


7 



+ u 


B DREG 


OJ 





(48) 
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and in DRED-scheme 


7 = 4 a 


(I + OM ) +a( ^ + o (1) ) +a2( - +ok 0) 


(49) 


Thus, in the framework of DRED scheme the vanishing contribution at j —> 1 is ob¬ 
tained for the NLO contribution to the anomalous dimension. 


3.6 Symmetry between 7 and j — 7 

Using the scale So = q 2 , the expression for u as a function of the correctly defined anoma¬ 
lous dimension 7 at general n can be written in the following form 


n 


^ = 4a( x (n ,7 — - ) + a 


5 ( n i 7 - J - X (n, 7 - ^J • X [n, 7 - | 


(50) 


where 


n 


X Ub 7 ~ 


5(n.7-f) 


2 4 /( 1 ) 4 /( 7 ) — 4/ (n + 1 — 7 ) ( 51 ) 

4/" ( 7 ) + 4 ’"(n + 1 - 7 ) - 2 $ (n, 7 - ^ — 2 $ fn, 1 — 7 + ^ 


6C(3) 


To calculate the anomalous dimension 7 we write oj in the ’’Lorentz invariant” form 


u = 4a ^2 4 /( 1 ) — 4/(7) — T(j — 7 ) + A( n, 7 ) <ij , j = n + 1 + a;, 


where 


A( n, 7 ) = S (n, 7 — ^ + 2 4 d( 7 ) + T'(n + 1 - 7 ) x (jb 7 ~ ^) • 

Using the analysis of the subsection 3.2, Eq. (p2p can be presented as follows 
u = 4a (2 4/(1) — 4 /( 7 ) — 4/(n + 1 + a; — 7 ) + e) , 
where £ can be written in the ” holomorphically separable” form 


(52) 


(53) 


(54) 


UJ 


£ = 


Here 


2 (p(7) +p{l + n- 7)) + n(0(7) + 0(1 + n - 7)), 
a; = 4a(2 4 r (l) — 4/(7) — 4 /(n + 1 — 7)) + 0(a 2 ) . 


1 00 1 1 

pin) = + *'(7) - oC( 2 ) = 2E 77 , o,w ~eC(2) 


( 55 ) 


(7 + 2 kf 2 


k=0 


and 0 ( 7 ) is given by Eq. 
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Let us calculate A(n, 7 ) near its singularities. To begin with, we consider 7 —> 0 for 
physical n > 0 : 

A(n, 7) -4 (^(!) - d»(n + 1 )) + -c(n ), 

7 7 

c(n) = 2 T'(n + 1) - 2 T'( 1 ) - (3\n + 1) - (3'{ 1). 

Therefore 

7 = ^ (l + u; ( T(l) - T(n + 1))) + ^ (U(l) - T(n + 1) + |c(n)) . 

At n = 0 the correction ~ a is absent, but for other n, especially for n —► —r — 1 , we have 
the large correction to 7 


A 7 = 4a (T(l) — T(n + 1)) , 

which has the pole at j —>■ —r and leads to a contribution changing even singularities of 
the Born term. It is related with the fact, that for positive n we calculate the anomalous 
dimensions of the higher twist operators (with the anti-symmetrization between n trans¬ 
verse and 1 + u longitudinal indices). 


Because in the case N = 4 SUSY the result is analytic in |n|, one can continue the 
anomalous dimensions to the negative values of \n\. ft gives a possibility to find the 
singular contributions of the anomalous dimensions of the twist -2 operators not only at 
j = 1 but also at other integer points j = 0, —1, —2.... In particular, as it was discussed 
above, in the Born approximation for the anomalous dimension of the supermultiplet of 
the twist-2 operators we obtain 7 = 4 a (T(l) — T(j — 1)) which coincides with the result 
of the direct calculations (see JT 8 |, [L9| and the discussion below). Thus, in the case N = 4 


the BFKL equation presumably contains the information sufficient for restoring the kernel 
of the DGLAP equation. 


4 The anomalous dimensions matrix in the = 4 
SUSY 

The DGLAP evolution equation for the moments of the parton distributions for N = 4 
SUSY has the form 

^l4g 2 A(j’Q 2 ) = XU ab(j)fb{j,Q 2 ) (a,b = q,g,(p), (56) 

A/a(i,Q 2 ) = XU ab(j)Afb(j,Q 2 ) (a,b = q,g), (57) 

where the anomalous dimension matrices 7a&(j) an 7 a b{j) can be written as expansions 
over the coupling constant a in the form: 

7 ab(j) = o- 7 a?(j) + « 2 - 7 SU), 7 ab(j) = a ■ 7^ (j) + ^ ’ 7 ^CO■ (58) 
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In the following subsections we will present the results of exact calculations for the 
leading order (LO) coefficients 7 ^ b \j) and 7 j$(j) and construct the anomalous dimensions 
of the multiplicatively renormalizable operators of the twist-2. In the NLO approximation 
the corresponding coefficients 7 ^(j) and 7 ^(j) are yet unknown and their calculation is 


in progress |26|. However, the form of the LO anomalous dimensions of the multiplicatively 
renormalizable operators in N — 4 SUSY is rather simple because they are expressed in 
terms of one function. Taking into account the universality of this result, the knowledge of 
the anomalous dimensions in the QCD case, the NLO corrections to the BFKL kernel jD|] 
and an experience in integrating some types of the Feynman diagrams, one can write an 
ansatz for the NLO anomalous dimensions of the multiplicatively renormalizable operators 
in the N — 4 SUSY. This result will be checked by us later with direct calculations of the 
matrix elements 7 ^(j) and 7 H\j). 


4.1 The results of exact calculations of the anomalous dimen¬ 
sions matrix in the N = 4 SUSY 

The elements of the LO anomalous dimension matrix in the N — 4 SUSY have the fol¬ 


lowing form (see ||T9|): 

for usual tensor twist -2 operators 

7 £>(i) = 4 0(1)- *0-1) 

2 2 



7 SHj) 


7 S(j) 


= 2 



TisC?) = 12 {jTl 


+ 1 J + 2 


7 = 


8 

j' 


7S?(i) 




4 («(i) - *(j) + j - jY 


Itpqij) = 


3 + 1’ 


7 $>U) = 4(T(1)-T(J + 1)), 
and for the pseudo-tensor operators: 


'40 i = 4 (Ay 


(59) 


= 4 *(i)-*y + i) 


7, 


S (0) W) = 



fJS’o) = 2 0- 


3 + 


7 ™(3) = 4 T(l)-T(j + 1) + - 


3 


(60) 


Note, that in N — 4 SUSY there are also the twist-2 operators with the fermion 
quantum numbers but the anomalous dimension of the corresponding multiplicatively 
renormalizable operators is the same as that for the bosonic components of the supermul- 
tiplet. 


13 


























4.2 Anomalous dimensions and twist-2 operators with a multi¬ 
plicative renormalization 


Let us denote the distributions of partons with the spin S by fs(x) and their correspond¬ 
ing momenta by fs(j)- We introduce also the inclusive probabilities A fs(x) and their 
momenta A fs(J) which are the differences of the distributions of the partons with the 
hclicities ±5. 

ft is possible to construct 5 independent twist-two operators with the multiplicative 
renormalization. The corresponding parton distributions and their LO anomalous dimen¬ 
sions have the form (see 0): 


fl U ) = fg(j) + fq (j ) + f<p(j)~ fq,g, V U ) > 

7 ?\ j ) = 4(T(1) — T(j — 1)) = —4S' 1 (j — 2) = 7+ ) (j), 

fllti) = -2(j-l)/ 9 (j)+/ g (j) + ^(j + l)A(j)~/gV^)’ 
ifhj) = 4 (^(l) — 4/(j + 1 )) = —4S' 1 (j) = 7o 0 ) (j), 

find) = - + fq(j) - ~ fq,g, V U)> 

7Ski) = 4 (T(l) — T(j + 3)) = —4S' 1 (j + 2) = 

= 2Af g (j) + Af q (j) ~ Af£ g (j), 

7 'iv(j) = 4 (T(l) - T(j)) = —451 (j - 1 ) = 7 +\j), 

fvU) = ~U~ tfAflti) + ~ A fq,gU)i 

7 fiU) = 4( v h(l) — T(j + 2)) = —4S' 1 (j + 1) = (61) 

Thus, we have one supermultiplet of operators with the same anomalous dimension 
7 LO (j), proportional to T(1) — T (j — 1 ). The momenta of the corresponding linear com¬ 
binations of the parton distributions can be obtained from the above momenta by an 
appropriate shift of their argument j in accordance to this universal anomalous dimen¬ 
sion 7 LO (j). Moreover, the coefficients in these linear combinations for N = 4 SUSY can 
be obtained from the conformal supersymmetry (cf. Ref [|17]]) and should be the same for 
all orders of the perturbation theory. 


4.3 The NLO anomalous dimensions and twist-two operators 
with the multiplicative renormalization 

We have the following initial information to be able to predict the NLO anomalous di¬ 
mensions of twist-two operators with the multiplicative renormalization in N=4 SUSY: 


14 






1 . As it was shown in the previous subsections LO anomalous dimensions are 
meromorphic functions. Moreover, there is really only one basic anomalous dimension 
7 LO (j) and all other anomalous dimensions can be obtained as 7 LO (j ± m), where m is 
integer. It is useful to choose: 


7 L °(j) = 4 (4/(1) - *(j - 1)) = -4 S^j - 2), (62) 

Then, 7 LO (j) has a pole at j —> 1 and vanishes at j — 2 . It is natural to keep the 
above universality also for the NLO anomalous dimensions 7 ^ (j) and 7 ^(j)- Adding 
some other properties we will construct ansatz for the basic NLO anomalous dimension 
l NLO (j). 

2 . There are well known results for NLO corrections to the QCD anomalous dimen¬ 
sions. 


3. In MS-scheme (with the coupling constant a) and also in the MS- like scheme with 
the coupling constant a (i.e. in the scheme based on DRED procedure), the terms rs./ C( 2 ) 
should be cancelled in the final result for the forward Compton scattering (see [p7fl - [p9| ). 
Therefore the terms ~ £(2) should be cancelled at even j in anomalous dimensions for 
the structure functions F 2 and F L (and for the unpolarized parton distributions) and at 
odd j in anomalous dimensions for structure functions g\ and F :i (and for the polarized 
parton distributions). 

4. From the BFKL equation in the framework of DRED scheme (see (|34|), (|k]) , ( |30| ) 
and m) we know, that there is no mixing among the functions of different ’transcen- 
dentality level’ [] i, i.e. all special functions at the NLO correction contain the sums of 
the terms ~ 1/n* (i = 3). Indeed, if we denote the following functions by the terms in 
the corresponding sums 

T ~ 1/n, 4T ~ f3' ~ C(2) ~ 1 /u 2 , h” ~ (3" ~ C(3) ~ 

then for the BFKL equation the LO term and NLO one have the ’levels’ i = 1 and i = 3, 
respectively. 

Because in N=4 SUSY there is a relation between BFKL and DGLAP equations, the 
similar properties should be valid for anomalous dimensions themselves, i.e. the basic 
functions 7 LO (j) and 7 NLO (j) should be of the types ~ l/j l with the levels i = 1 and 
i = 3, respectively. The LO basic anomalous dimension is given by Eq. (^). Then, the 
NLO basic anomalous dimension r y NLO (j ) can be expressed through the functions: 

Si(j ~ 2), Kj(j - 2 ),S ki i(j - 2 ),K kil (j - 2) : gk)S l (j - 2 ), ((k)K i (j^2), ((i) 

4 Note that similar arguments have been used also in p0| for finding results for some types of compli¬ 
cated massive Feynman diagrams. Their evaluation was based on the direct calculation of several terms 
in the series over the inversed mass and on the knowledge of the basic structure of the series obtained 
earlier in |[iTT| . |yi|i by an exact calculation of few special diagrams with the use of the differential equation 
method fl32fl. 
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3 and j + l — 

- i), where 

SiU) = 

t Y ~ t'-'o + 1). 

S2,i(J) = 

,7Zi m 

K 2 (j) = 

(1 - (-iy)|c(2) + <- 

KsU) = 

(1 - (_iy)| C (3) + (- 

Agi (j) = 

(1 - (-iy)|c(3> + <- 


m= 1 
3 


m— 1 


m=l 


m 


(-i) 


m +1 
2 

m+1 




m 


P (j + 1 ); 

/(j + 1), 


(63) 


Note that the terms ~ £(2) should be absent in accordance with item 4. 
Moreover, the terms 

S,U-2)/(j±m)\ Ki(j — 2)(j ± m) k (j + l = i) 


(64) 


should be absent, too. There are two reasons for this absence. 

First one, these terms have additional poles at different values of j: j = =p m. But 
such additional poles should be absent, if we start with the BFKL equation and obtain 
l NLO {j) by replacement n —► — 1 — r because the BFKL equation does not dublicate the 
poles at different values of j. 

The second reason comes from the consideration of linear combinations (|6T|) . If, for 
example, in the polarized case terms (Q) contribute, then we will have the terms ~ 
(j ± m)^ 1 in one combination and the terms ~ (j ± m ± 2) -1 in another combination. 
However, from the direct calculation of QCD NLO anomalous dimensions in the polarized 
case (see |33]]) we know that only the terms ~ j^ 1 and ~ (j +1) _1 can contribute to these 
combinations. 

So, terms (^4j) should be absent in the results for the universal NLO anomalous di¬ 
mension in the N=4 SUSY case. 


5. Further, the NLO anomalous dimension y ArLO ( j ) is equal to a combination of 
more complicated contributions (i.e. the contributions containing the functions with the 
maximal value of i: i = 3) for the QCD anomalous dimensions (with the SUSY relation 
for the QCD color factors Cf = Ca = N c ). 

Note, that these most complicated contributions (with i = 3) are the same for all 
QCD anomalous dimensions with above SUSY relation (with a possible exception for the 
NLO scalar-scalar anomalous dimension which is not known yet). 

Thus, for N = 4 SUSY the NLO universal anomalous dimension 7 NLO (j) has the form 

l NL °(j) = 16 Q(j — 2), (65) 
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where 


QU) = K a , l (j) + ^S 3 (j)-K s (j))+S 1 (j)(s 2 (j)-^a(j)) (66) 

We would like to note that the sums K 2 (j), K 3 (j ) and K 2 .i(j) (see Eq.(|63D) have been 
calculated directly only at even values of j (or at odd ones, when (—1) J —>► (—1) J+1 ). The 
analytical continuation to complex values of j can be done rather easily (see |25|] and 
respectively). 

6. We can add the term ~ £(3) to the r.h.s. of (j 66 |), but due to the condition 
7 NLO (j = 2 ) = 0 it cancels. 

So, for N=4 SUSY the universal anomalous dimension 7 (j) has the form 

7(i) = S7 LO (j) + a 2 'y NLO (j), (67) 

where 7 LO (j) and 7 NLO (j) are given by Eqs. ( |62D and (|65|), respectively. All other anoma¬ 
lous dimensions can be obtained as 7 LO (j ± m) and 7 NLO (j ± m), where m is an integer 
number. 



Thus, the above arguments allow us to construct the NLO corrections to anomalous 
dimensions in the N=4 SUSY, which were unknown earlier. We plan, however, to check 
these results by direct calculations [^ 6 j. 


4.4 DGLAP evolution 

Using our knowledge of the anomalous dimensions we can construct the solution of the 
DGLAP equation in the Mcllin moment space in the framework of N=4 SUSY. 

A. Polarized case 

The polarized parton distributions are splitted in the two contributions: 


= A/+(j,Q 2 ) + A/-(j,cn 


( 68 ) 


where 
at LO 


( 0 ) 

Q 2 ) = A/JT°(.7 Ql ) ('7± ) = 4Si(j =F 1)) 


(69) 


at NLO 


A 4 ± 3 Cb Q 2 ) = A f£g NLO U,Q 2 0 )(^) 


2 \ (7± 0) o+7±±a A 


( 7 E = 16 QU T 1)) , (70) 
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where 


Mt; 9 NLO 


U, Q o) 


+ 


A f±,NLO 
J q,g 


(i, Ql) 


7±i a \ 

#-7 


7Si a 


~ (0) 

7V 


~(o) 
- 7 : 




(J, Ql) 


( 71 ) 


Notice that only the anomalous dimensions are important for N = 4 SUSY at the 
order 0 (a 2 ): they contribute to the (Revolution of parton distributions. 

The anomalous dimensions y£y give contributions at 0{a 2 ) only to the normalization 
factors A f^g NLO (j, Ql)- Their contribution to the (Rdependent part of Af± g (j, Q 2 ) starts 
at 0 (a 3 ) level in the following form: 


U’Ui 


~ (o) 
7± 


~ (o) 
7t 


(72) 


B. Nonpolarized case 

The polarized parton distributions are splitted in the three parts: 


f m u,Q 2 ) = E /b«.(i,e 2 ), 

,0 


(73) 


whereQ 
at LO 


fLMQ 2 ) = 


(o) 

2 \ Ti a 


(t f > =4S,(j-i2)) (74) 


at NLO 


r„,jj,Q 2 ) = r,Z°u,Ql)[^ 


2 \ 


( 7 ? ) = 16gC/-i2)) (75) 


As in the previous case A, only anomalous dimensions yj^ are important at 0(a 2 ) 
in N=4 SUSY. The anomalous dimensions 7 ^ (i 7 ^ k) contribute at 0(a 2 ) level only to 
normalization factors fg^°(j,Q o)- 


5 Conclusion 

Above we reviewed the results for the next-to-leading corrections to the kernel of the 
BFKL equation and to the anomalous dimensions of twist-2 operators in the extended 

5 The arguments are j — i2 = {j — 2. j. j + 2} for respectively i = {+,0, —}. 
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N — 4 SUSY. The absence of the coupling constant renormalization in this model leads 
presumable to the Mobius invariance of the BFKL equation in higher orders of the per¬ 
turbation theory. The cancellation of non-analytic contributions proportional to <5° and 
in N — 4 SUSY is remarkable (such terms contribute to u in the framework of QCD 
as it has been demonstrated in |J |35|). This property could be a possible manifestation 
of the integrability of the reggeon dynamics in the Maldacena model |3B| corresponding 
to the N — 4 SUSY in the limit N c —> oo. Note, that in this model the eigenvalues of the 
LLA pair kernels in the evolution equations for the matrix elements of the quasi-partonic 
operators are proportional to — 1) — ^(1) [|l^, [19|, which means, that the correspond¬ 


ing Hamiltonian coincides with the local Hamiltonian for an integrable Heisenberg spin 
model. The residues of these eigenvalues at the points j = —k are obtained from the 
BFKL equation by an analytic continuation of the anomalous dimensions to negative in¬ 
teger values of the conformal spin \n\. Therefore the DGLAP equation is not independent 
from the BFKL equation in N = 4 SUSY and their integrability properties at N c —> oo 
are presumably related. 
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